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The impact of radio-frequency (rf) currents on the direct current (dc) driven switching dynamics in 
current-perpendicular-to-plane nanoscale spin valves is demonstrated. The rf currents dramatically 
alter the dc driven free layer magnetization reversal dynamics as well as the dc switching level. This 
occurs when the frequency of the rf current is tuned to a frequency range around the dc driven 
magnetization precession frequencies. For these frequencies, interactions between the dc driven 
precession and the injected rf induce frequency locking and frequency pulling effects that lead to a 
measurable dependence of the critical switching current on the frequency of the injected rf. Based 
on macrospin simulations, including dc as well as rf spin torque currents, we explain the origin of 
the observed effects. 



INTRODUCTION 

Direct electrical currents, via the spin transfer torque 
(STT), can excite a broad variety of magnetization dy- 
namics, including reversal and steady-state precessional 
states in current-perpendicular-to-plane spin valves. 
These precessional dynamics, stabilized by the oppos- 
ing effects of the STT and the damping torque, in ad- 
dition to being of great scientific interest have important 
implications for the development of new radiofrequency 
(rf) devices Q, i]. While prior studies focused on the 
dynamics in the non-hysteretic regime, it has also been 
recognized that precessional dynamics can exist in the 
hysteretic regime @, 0, S 0- These pre-switching (PS) 
precessional modes are of great interest as they are re- 
sponsible for driving the switching dynamics and control 
of switching is key to potential memory and logic appli- 
cations. 

In this work we investigate the highly unexplored PS 
dynamical regime as well as the effects of additional rf 
currents on the PS and switching modes. The PS modes 
are found to interact with the applied rf currents having 
frequencies in the range of the dc-only driven modes, and 
under certain conditions are found to frequency lock. In 
contrast to previous frequency locking results @, @, H|, 
by looking at these effects in the PS regime, we find a 
measurable dependence of the dc switching threshold on 
the external rf frequency, fi n j . In the locking range the rf 
appears to have a stabilizing effect which hinders switch- 
ing, however, at frequencies just below this range, we find 
a measurable reduction in the direct critical current for 
switching, I c . While we have already reported this effect 
[l(]| , an explanation for the origin of this effect has not 
been given. In this work, by studying the behavior of this 
system in frequency domain and through macrospin sim- 
ulations (MS) (using the Landau-Lifshitz-Gilbert equa- 
tion with Slonczewski STT 0]) we provide a description 
of the magnetization dynamics originating from the in- 
teraction between the dc driven switching modes and the 
applied rf. The results show good agreement with our 
experimental data. 



EXPERIMENT 

The samples are 50 nm x 100 nm hexagonal pillars 
comprised of an IrMn pinned antiparallel (AP) coupled 
bilayer Co 50 Fe 5 o(25A)/Ru(8A)/Co 50 Fe 5 o(25A), a 40 1 
Cu spacer, and 35 A Ni 92 Fes as the free layer (FL). 
We assume all observed current induced magnetic ex- 
citations occur in the FL. In our system of reference, 
positive fields/negative currents stabilize the AP mag- 
netization state. Here we present results from a single 
element, selected due to its comparatively high power 
PS modes even in absence of applied external field H app . 
This is likely due to the significant dipolar field Hdi P ~ 
170 Oe, at the free layer (from the field hysteresis curve 
offset), despite the use of AP coupled pinning layers. We 
point out, however, that variations from one element to 
another are observed and furthermore, PS dynamics are 
not detectable for all pillars. Nevertheless, qualitatively 
similar results were obtained for all samples for which 
similar amplitude PS modes were observed. 

Experiments were performed in a variable temperature 
Desert Cryogenics probe station equipped with a super- 
conducting magnet for applying in-plane magnetic fields. 
Contact to the pillar was made through high frequency 
(<40GHz) probes connected to the (rf+dc) mixing port 
of a broadband bias tee. The inductive port of the bias 
tee was used for low frequency dynamic resistance mea- 
surements using a lock-in amplifier, while the capacitive 
side was attached to a splitter that reroutes rf signals 
from a high frequency signal generator and toward a + 
44 dB amplifier attached to a spectrum analyzer. While 
the sample stage was held at 4.2 K, the true sample tem- 
perature, when injecting direct currents, I dc on the or- 
der of 1 mA, was estimated to be approximately 25 K 
due to Joule heating. All experiments were performed 
at low temperatures in order to minimize thermal fluc- 
tuations. Calibration of the electrical set-up was carried 
out by measuring the temperature dependent Johnson 
noise from 50 Ohm test resistors, and a network analyzer 
was used to measure the complex and frequency depen- 
dent impedance Z$ of the sample, which was then used 
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to estimate the rms amplitude of the injected rf current. 



RESULTS AND DISCUSSION 

The macrospin model predicts the onset of sub-critical 
current driven stable precessional dynamics around the 
local dipolar field for currents above a threshold value I t , 
for which the STT balances the damping torque. As I dc 
increases, the orbits are predicted to expand while pre- 
cession frequencies decrease 0, @. Experimentally, we 
identify potential regions of STT induced dynamic insta- 
bilities in the FL magnetization as anomalies in dV/dl. 
Fig. 1(a) illustrates this for current sweeps taken at H app 
=25, 50 and 100 Oe, where the field is applied a few de- 
grees off the easy axis direction. We concentrate mainly 
on the STT magnetization dynamics prior to the AP to P 
transition, occurring approximately at -560 Oe at I dc =Q 



(not shown) and 1.45 mA at H, 



app~ 



=0. 



At currents greater than the instability threshold cur- 
rent It, anomalies, appearing as small dips for H app =25, 
50 Oe (hysteretic) and larger dips for H app —100 Oe (non- 
hysteretic) in Fig. 1(a), reveal strong microwave activity. 
Shown in Fig. 1(b) is a phase diagram, where the color 
scale represents dV/dl for I dc swept from negative to 
positive values and the continuous lines show the switch- 
ing boundaries obtained for current sweeps in both direc- 
tions. The PS regime appears as a narrow bright area, 
on the right side of the dotted line and to the left of the P 
region. Low positive fields (<100 Oe), that oppose STT 
yet sustain hysteretic behavior, stabilize and strengthen 
these precessional states over a wider I dc range, as can 
be seen through the widening of the PS regime. Fig. 
1(c) compares current driven dV/dl loops measured with 
and without injected rf at H app = 0, showing the strong 
anomalies induced by the rf. We have observed this type 
of behavior for samples with two distinct FL materials 

Examples of power spectral density (PSD) data in the 
PS regime are shown in Fig. 2(a). As I dc increases in this 
range, the emitted power rises as the dc driven resonance 
frequency /q c drops (redshifts) in agreement with the 
macrospin model. In Fig. 2(b) we plot the measured /g c 
vs I dc for H app —0 and 50 Oe, showing the peak frequen- 
cies decrease as the current rises toward the switching 
value I c . Extrapolation of these data through I c suggests 
that the minimum resonance frequencies approaching re- 
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and for H avv —0 are near 4.2 GHz. Fol- 



lowing Kiselev et al. [3j , we make a macrospin estimation 
of the precession and misalignment angles (between the 
pinned layer and the FL) by assuming small angle preces- 
sion and measuring the total integrated power in the first 
and second harmonic. While this approximation may be 
imprecise, in view of the highly non-Lorenztian spectra, 
it provides a rough picture of the behavior of the pre- 
cessional dynamics as I dc is ramped. We find precession 
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FIG. 2: (Color online) (a) Power spectral density measured at 
H apP =50 for I dc values below the switching value. The PSD 
measured with I dc =0 has been subtracted, (b) /q c vs. I dc 
for H app =0 (squares), 50 Oe (circles) (symbols at peak value, 
bar denotes FWHM span). 



angles increasing from 10 to 30 degrees, for I dc increas- 
ing from 1 to 1.4 mA and for < H app < 100 Oe. How- 
ever, the multiple peak structure and broad lincwidths, 
reaching up to ~ 500 MHz, suggest incoherence and co- 
existence of several modes. Such incoherent behavior has 
been observed in micromagnetic modeling [Tlj |. an d can 
be understood from the spatially non-uniform local de- 
magnetizing fields at the edges of the spin- valve. 

As previously stated, the effect of adding an rf current, 
I r > to I dc on the FL STT dynamics can be seen from the 
dV/dl vs current hysteresis loops in Fig. 1(c). This fig- 
ure includes data for the no rf case (dashed line) and for 
three different injected frequencies fi n j =4.3, 4.7 and 5.5 
GHz. In all cases, a power of -15 dBm was applied at the 
signal generator, resulting in (/ r /) rms ^ 1mA. The dips 
(or peaks) appearing before the AP to P (or P to AP) 
transition are caused by interactions between and 
the I dc driven precession resulting in changes in the FL 
magnetization precessional trajectory and in the magne- 
toresistance. Focusing on the AP to P switch, note that 
as finj increases, the I dc level at which the dip occurs de- 
creases, following the /q c increase induced by decreasing 
I dc . As discussed below, for external frequencies below 
the /g c range, we observe a reduction in I c , 
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FIG. 3: (Color online) (a) Comparison of the PSD with (cir- 
cles, squares) and without (star) rf for I dc in the PS range 
and Happ—25 Oe, showing frequency locking for finj —4.5 
GHz and 4.75 GHz. (b) R dc for I dc near the AP to P switch, 
and for several values of injected frequencies. The dashed line 
shows the no rf case. 

To elucidate the effects of the rf on the FL dynamics, 
we studied the changes in the frequency domain behav- 
ior induced by adding I rf to I dc . Results suggest partial 
frequency locking between the dc driven precession and 
the external rf signal as shown in Fig. 3(a). For these 
experiments power was limited to -23 dBm to avoid sat- 
urating the amplifier. Low applied fields below 50 Oe 
were applied to enhance the signal originating at the sam- 
ple. Fig. 3(a) illustrates the partial frequency locking of 
the applied and internal oscillations by comparing signals 
obtained at 7 dc =1.40 and 1.44 mA, with (after subtract- 
ing the external rf signal) and without external rf. While 
the peak frequency for l r f=0 is at 4.8 GHz, applying 
with fi n j at 4.5 GHz or 4.75 GHz narrows the signal and 
shifts the measured power to the fi n j value. As expected, 
locking only occurs when I r f is applied at a frequency 
at which the dc driven mode also has significant power. 
Further measurements (not shown here) suggest a corre- 
lation between the existence of dips/peaks in dV/dl and 
the (I dc , finj) values for which frequency locking occurs. 

Additional insight can be extracted from dc resistance 
(R dc ) measurements. As shown in Fig. 3(b), for fi n j close 
to /g c , a reversible plateau forms prior to switching, with 
a AR of approximately 0.05 f2. The effect appears at 
lower values of dc for higher injected frequencies. Our 
data confirms that the (I dc , finj) values for which the 
dips/peaks in dV/dl form and for which frequency lock- 
ing occurs also coincide with the formation of the plateau 
in the dc resistance. In accordance with the previously 
observed behavior, the Rdc data shows that for fi n j in 
the range (J^ c ) mm < 4.2GHz, the plateau does not form 
and I c drops as the fi„j increases. 
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FIG. 4: (a)(J dc , f inj ) phase map obtained experimentally at 
H apP = with R dc represented through gray scale. The white 
squares show the locations of the measured /q c . (b) (I dc , finj) 
phase map obtained through MS, with gray scale representing 
average final state MR. The white squares show the locations 
of the /o c obtained numerically. 

In Fig. 4(a) we plot R dc in gray scale for values near 
the AP to P transition, thereby generating an (I dc , finj) 
phase diagram. The diagram shows the reduction trend 
in I c as fi 7l j increases in the range below (/o c ) mm ~ 4.2 
GHz. This reduction in I c is well above the sweep-to- 
sweep variation. For this specific sample, we estimate 
the magnitude of the variations in I c by measuring ~ 80 
loops with finj— 4 and 5 GHz. We account for a slight 
frequency dependence of the sample impedance by apply- 
ing at each fi n j , the power required to pass an equal mag- 
nitude (j r fy ms = i m A through the sample. The differ- 
ence between the average switching currents < if GHz > 
and < GHz > j g g rea ter than 15%. The standard de- 
viations on the other hand, are below 3%. 

For fi n j in the /q c range, we identify the plateau, la- 
beled as the frequency locking region, as well as an in- 
crease in I c . The open squares show the approximate 
locations of the f dc (I dc ) peak resonances, which are cen- 
tered in the locking region as expected. While the ap- 
pearance of the plateau in R dc could be attributed to the 
MR effect resulting from an rf induced modification of 
the dynamical trajectory of the magnetization, frequency 
locking between the external rf current and the dc driven 
precession is also expected to produce a dc rectification 
voltage 0, [H . Rectification effects due to frequency lock- 
ing have been studied in similar samples, when appearing 
in the non-hysteretic regime 0,H, 3- 

Further physical insight into the rf driven effects can be 
gained through MS based on the Landau-Lifshitz-Gilbert 
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equation using Slonzcewski @ type STT. We compare 
simulation results for the cases of T rf =0 and I rf ^ 0. 
For an initial AP configuration and no rf, and using 
a=0.025, we observe switching for I dc > 1.4 mA within 
a 10 ns window. For lower I dc we find PS dynamics. The 
calculated Fourier spectral density of the in-planc com- 
ponent of the magnetization along the short axis, M y (t) 
drop from 4.1 GHz to 3.1 GHz as the current increases 
between 1.2 to 1.3 mA. The numerical experiments also 
provide information concerning the effects of I r f+ dc on 
the dynamics in time and frequency domain. An im- 
portant conclusion we extract from the MS is that the 
precessional frequency and orbit are correlated such that 
small (large) orbits correspond to high (low) frequencies 
(which is also the case for dc-only driven PS precession). 
Therefore, when I r f with fi n j ~ /q c is applied such that 
frequency pulling or frequency locking occur, effects on 
the magnetization precessional orbits also appear. As 
the MS do not include the Oersted field from I dc + r f ; we 
conclude that all rf induced changes on the magnetiza- 
tion trajectory originate from STT. However, rf magnetic 
field induced reductions in switching fields have also been 
observed [11 El- 

In Fig. 4(b), we show the T = simulated switching 
boundary by plotting the MR (average value for the fi- 
nal 5 ns of a 20 ns simulation window, normalized to the 
experimentally obtained MR amplitude) as a function of 
I dc and f lnj . The rf current was set to (I r f) rms =1 m A 
to match the experimental amplitude. The light/dark 
areas correspond to high/low MR states. The switch- 
ing boundary shows a reduction in I c for fi n j ~ 3 GHz, 
which is close to the numerically obtained /q c , for I dc 
near I c . At slightly higher fi n j ~ 3.5 GHz, I c reaches 
a maximum. A slow decay follows as fi„j continues to 
increase above 3.5 GHz. The gray area labeled C shows 
a drop in the MR with respect to the AP state. When 
passing from region B to region C, by increasing fi n j at 
a fixed I dc , the simulated PS orbits show an an abrupt 
expansion above 10°. This expansion appears as a sharp 
contrast between regions B and C, revealing the bound- 
ary for a major change in the induced precessional orbits. 
The Fourier transform spectrum of M y shows that the 
dynamics within the gray area (C) contains a dominant 
spectral component at finj- The scattered gray points 
appearing at I dc > I dc for f- m j ~ 3 GHz, correspond 
to sustained high amplitude MR oscillations. Additional 
MS (not shown) reveal that these gray points, together 
with other gray areas near the switching boundary, be- 
come dark (switch) at higher temperatures, thereby ac- 
centuating the I c minima occurring between 3 and 3.5 
GHz. Furthermore, the scattered points observed at T=0 
may be a numerical artifact, as their locations depend on 
the initial conditions. Hence, the reduction for /j n j below 
/q c and the increase for f in j ~ /q c , in 7 C , as well as the 
frequency locking boundary are qualitatively reproduced 
by the simulations. 



To understand how the rf affects the switching bound- 
ary, recall that as the I dc is ramped toward the switch- 
ing value, the corresponding resonance frequency drops 
toward (/o c ) mm - If the injected current additionally 
has an rf component with frequency fi n j just below 
(/g c ) mm , then frequency pulling will open (f in j pulls to- 
ward a lower frequency) the orbit as I dc approaches I c . 
This pull drives the system toward finj, however, since 
finj < (/o c ) mm , there is no stable orbit corresponding to 
finj , thus frequency locking does not occur and switching 
at a reduced I c is observed. This effect occurs at lower 
I dc values for higher fi n j in this range and this may be 
the origin of the consistent decrease in I c obtained as fi„j 
increases in the range below (/p c ) mm - 

In contrast, for f in j > (fo c ) mm , frequency locking oc- 
curs for I dc near the switching boundary, stabilizing the 
orbit and increasing I c . As long as the precessional fre- 
quency is locked to f in j , the orbit does not expand as I dc 
increases. Therefore, switching does not occur until I dc is 
sufficiently high and the corresponding f d c is far enough 
from finj for the system to unlock and the magnetization 
to switch. For higher frequencies, such that f in j does not 
remain locked as I dc approaches J c , there are no expected 
effects (of this type) of the rf on I c , as is the case for even 
higher injected frequencies fi n j > (f dc ) max for which no 
frequency locking occurs. 

In conclusion, we demonstrate that rf currents tuned 
to the frequency range of the dc driven magnetization 
precession frequency range, can significantly alter the dc 
driven PS and switching modes. Based on frequency do- 
main measurements and MS that well describe our ex- 
perimental data, we explain the origin of the rf driven 
effects on the precessional orbits and frequencies and on 
the dc critical switching boundary. This mechanism for 
reducing I c by additionally injecting rf currents with fre- 
quencies slightly below (/Q C ) mm may be of technological 
interest for the design and control of STT devices. 
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